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Abstract 

We review the stocliastic Gross-Pitaevskii approach for non-equilibrium finite temperature Bose gases, 
focussing on the formulation of Stoof; this method provides a unified description of condensed and thermal 
atoms, and can thus describe the physics of the critical fluctuation regime. We discuss simplifications 
of the full theory, which facilitate straightforward numerical implementation, and how the results of 
such stochastic simulations can be interpreted, including the procedure for extracting phase-coherent 
('condensate') and density-coherent ('quasi-condensate') fractions. The power of this methodology is 
demonstrated by successful ab initio modelling of several recent atom chip experiments, with the important 
information contained in each individual realisation highlighted by analysing dark soliton decay within a 
phase-fluctuating condensate. 
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I. INTRODUCTION 



Many theories have been devised to study the static and dynamic properties of weakly inter- 
acting, ultracold, atomic Bose gases at finite temperatures 111121. An important feature of these 
systems is that a condensate coexists with a non-condensed component beneath the temperature 
for the onset of Bose-Einstein condensation (BEC): repulsive interatomic interactions cause a 
depletion of atoms from the condensate, while thermal effects additionally promote atoms from 
the ground state of the system. An accurate description of partially condensed Bose gases thus 
requires a theory capable of describing condensed and non-condensed fractions within a unified 
framework. Symmetry breaking finite-temperature approaches rely on the explicit existence of a 
condensate mean field; such approaches are therefore useful away from the critical region, where 
there is a well-established condensate. Nonetheless, at temperatures close to the transition point, 
critical fluctuations make the definition of a condensate mean field much more difficult, an impor- 
tant consideration when studying condensate growth, or low dimensional geometries, for which 
the temperature range over which fluctuations are important is broader. Fluctuations play a key 
role in such cases and their presence motivates a stochastic description of the weakly interacting 
Bose gas. While several such methods are discussed in this book, this Chapter focuses on the 
non-equilibrium formulation of Stoof [!3]-(7l, whose end result is a nonlinear Langevin equation for 
the field representing the condensate and its fluctuations. 

II. METHODOLOGY 

To treat the fluctuations inherent to the process of Bose-Einstein condensation, a non-equilibrium 
probabilistic theory is required. Fokker-Planck equations, originally introduced to describe the 
Brownian motion of particles, achieve this by describing the time-evolution of an appropriate prob- 
ability distribution. The mapping of a Fokker-Planck equation to an equivalent representation in 
terms of a Langevin equation, is often made for numerical impementation; such equations appear 
frequently in diverse fields including financial market modelling [[8]|, superconductors ||9l [TOl . 
high-energy physics [fTTll and turbulence [12J. 
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A. Fokker-Planck Equation for an Ultracold Bose Gas 



In Ref. jH, Stoof derived an equation for the evolution of the full probability distribution for 
the weakly-interacting Bose gas using a field-theoretic formulation of the non-equilibrium Keldysh 
theory [[T3]| . within the many-body T-matrix approximation. Using a Hartree-Fock-like ansatz, the 
total probability distribution was separated into a product of respective probability distributions 
for the condensate and thermal particles, which led to two coupled equations for these 'subsys- 
tems'. Firstly, by integrating out the thermal degrees of freedom, the dynamics of the condensate 
distribution function, /'[O*, O; t], was found to obey: 



ot 



- J Jr ^^|-^V2 + yext(r)-/m0 + ^|O(r,0p-M0jo(r)P[O*,O;?] 
^ / ^r^^(-^V2 + yext(r) + /m0+^|O(r,0p-M0jo*(r)P[O*,O;?] 



2 J ~" (50(r)(50*(r) 



where, represents a functional derivative with respect to the complex field O. The term 

zi?(r, t) describes gain or loss of particles due to collisions which transfer atoms between the con- 
densate and thermal cloud given by 

X 6(S,. + ~S2-h- Sa) [flih + l)(/4 + 1) - (/2 + miA] , (2) 

where e,- = |p,p/2m -I- Vext(r) + 2^(|0(r, Op) is the Hartree-Fock energy of a thermal atom and 
fi = f(si, t) the Wigner distribution function for thermal atoms. 

The strength of the fluctuations is set by the Keldysh self-energy, defined by 

" = -^-■^^ / / (Is? / - - 

X 6{S, + ~S2-h- §4) [/2(/3 + l)(/4 + 1) + (/2 + DMa] • 0) 

The R term arises as the difference between the rate of scattering processes into and out of the 
low-lying modes of the system, R = n(r°"' - P")/2, where r°"' oc f2(fj + l)(/4 + 1) and P" oc 
ifi + l)/3/4> while the self-energy is the sum of these rates = -/^(F™' -t-P"). Thus, although at 



equilibrium scattering processes should be zero on average (r°"' ^ P"), fluctuations nonetheless 
persist, highlighting the dynamical description of the equilibrium state. 

The above Fokker-Planck equation for the condensate is coupled to a quantum Boltzmann 
equation for the distribution function describing the thermal cloud, 

^ + (Vpg) • (V/) - (Ve) • (Vp/) = Cnif] + C22U], (4) 

arising from the corresponding probability distribution evolution for thermal atoms. Scattering 
processes which transfer atoms between the condensate and thermal cloud are described by the 
coUisional integral 

X + €2 -S3- sdi^ntif [6(-p - P2) - 5(p - P3) - 5(P - P4)] 

X[(/2 + l)/3/4-/2(/3 + l)(/4+l)], (5) 

while thermal-thermal collisions are represented by 

^''^^^ = ^J (2^'J(2m^J (W^^^^^ ^(P + P2-P3-P4) 

X 6(8 +~S2-h- SaW + l)(/2 + l)/3/4 - //2(/3 + l)(/4 + D]- (6) 

Thus, this formalism may be interpreted as a stochastic number-conserving generalisation of the 
ZNG kinetic theory [fT4l[T5l . which supplements dissipative processes affecting the condensate, 
by essential fluctuations. 



B. Formulation as a Langevin Equation 

The energy Sc which appears in the expression for the self-energy and the damping term iR{r, f) 
is the energy associated with removing an atom from the condensate, which is given by the oper- 
ator (IMll 

= -;?-V' + VextW + ^|0(r, tf. (7) 
2m 

As discussed by Duine and Stoof in Ref. [jTil, the fact that is an operator dependent upon O leads 
to a complicated stochastic equation with multiplicative noise. To proceed, we may approximate 
the Wigner functions, /(r,p, 0, by Bose-Einstein distributions, henceforth denoted by nBE(r, p)- 
This assumption is consistent with an equilibrium thermal cloud, which we thus represent as a heat 
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bath with a chemical potential, jj., and temperature, T . This process leads to a relatively simple 
expression linking the damping term and the Keldysh self-energy [l5l|7]|, 

- iR{r, t) = i^2^(r, t) [hbeCsc) + 1/2]"' . (8) 

A key observation is that this represents the fluctuation-dissipation relation for the system: it 
describes the relationship between the magnitude of fluctuations, set by ^Z^(r, t), and the damping 
due to the source term zi?(r, t). This relation depends upon the equilibrium mode populations, set 
by the sum of thermal Bose-Einstein populations hbe, and an extra quantum contribution of half 
particle per mode, on average; physically, these respectively represent stimulated and spontaneous 
contributions to the scattering rate. As high energy atoms are here assumed to be close to thermal 
equilibrium, Eq. ([8]) should be valid in the regime of linear response, applicable to perturbations 
which do not strongly aff"ect the thermal cloud. 

The Fokker-Planck equation for the condensate, Eq. ([!]), can be mapped to an equivalent repre- 
sentation as a Langevin equation (see flU), known in this context as the Stochastic Gross-Pitaevskii 
Equation (SGPE), which takes the form 



..5O(r,0 
in- 



+ Kxt(r) - iRir, t) + gmr, tf - ^ 

2m 



O(r,0 + /7(r,0. (9) 



dt 

This can be identified as a T > generalisation to the usual Gross-Pitaevskii Equation (GPE); 
it includes the scattering of particles into/out of the thermal cloud (-/i?0), with condensate fluc- 
tuations modelled via the dynamical noise term tj. In order to obtain an equation that can be 
easily solved numerically, and upon noting that /3{Sc - yu) is small at high T (J3 = l/k^T), or close 
to equilibrium (Sc « //), we Taylor expand the Bose-Einstein distribution of Eq. ([8]) in terms of 
this variable. The result of retaining the leading order term in this expansion is to replace the 
fluctuation-dissipation relation of Eq. ([8]) with its classical counterpart, based upon the Rayleigh- 
Jeans distribution, yielding [l5]-(7]l 

- iR(r, t) = ^ni.\r, t) (£c - //) . (10) 

Using Eq. Q and Eq. ^ in Eq.Q we finally obtain 
. 5O(r,0 / 13 



_^V2 + yext(r)+g|O(r,0l'-// 
2m 



O(r,0 + /7(r,0, (11) 



with Gaussian noise ensemble correlations {ri*{Y, t)rj{r', t')) = /(^^/2)S^(r, t)6{r - r')5{t - 1'). This 
is the form of the SGPE solved to date [|6l [T6] - l24l . By analogy to other stochastic methods, 0(r, t) 



should now be understood as representing a unified description of atoms within the low-energy 
modes of the gas, up to some energy cutoff, that are in contact with a heat bath made up of the 
remaining higher energy thermal atoms. 



C. Stochastic Hydrodynamics 

The corresponding stochastic hydrodynamic theory led to the following generalisation of the 
usual continuity and Josephson equations |7|, 



^Hfpll + V . (n,(r, 0vs(r, t)) = - ^ill^ir, t) iiu,(r, t) - /u) n,(r, t) 
ot 2 

+ 2 ^jn,{r, t)^{r, t). 



(12) 



T«^2 (r, t) -— =11- iu,(r, t) + , (13) 

respectively, where and 6 are the condensate density and phase, judr, t) = -(fi^V^ ^Indv, t)/2m ^fnjr^)+ 
Kxt(r, t) + gndr, t) + (m/2)v^(r, t) and the velocity is given by the gradient of the phase Vs(r, t) = 
hVOir, t)/m. The noise terms are now given by <v(r, t)v(r', t')) = i{fp- IA)fiL'^{r, t)5{r - r')6{t - t') 
and {^{r,t)^{Y' ,t')) = i(L^{r,t)IA)6{r - r')6{t - t'). This approach is easily amenable to analytic 
variational methods fTTilSSl. 



D. Simple Numerical Implementation 



The numerical solution of the SGPE, Eq. ( [TT] ), is not much more complicated than the usual 
GPE: following Bijlsma and Stoof [l6ll26l, we seek to propagate the system from a time by a 
time step A?, via 



0(r, t„ + At) = exp (-/[I + /3ni:''(r)/4][e, - iu]At/n) 0(r, - (^)^m(r) 



(14) 



where we have defined the noisy field at the time step as ^m(r) = exp{-/[l + fiKE,^ (r) / 4]{Sc - 
li)tm/n}l^fl"'^^' df exp{/[l +j3ni.^(r)/4-][£,-pi]t'/n}^T](r,t'). This has correlations given by <^;(r)^„(r')) 
i{fi^/2)'Z'^{r)6ir - r')d,„nAt + 0(Af). Making use of Cayley's form for the exponential of Eq. ([14]), 
the problem is then reduced to the solution of 



|l + / [l +;0^2^(r)/4] (ec,m+i/2 - yu) At/2n} 0(r, t,,, + At) 

= [\- i[\ +;e^S^(r)/4] (e,,„+i/2 -//)Af/2;i) 



0(r,f™)-(-)^,„(r) 



(15) 
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where £c,m+i/2 = [sdtm) + ec(tm+\)] /2 is the operator of Eq. (|7]) evaluated at the mid-point of the 
time step [|26l : the remaining step is to spatially discretise Eq. ( [T5] ), which can be achieved using 
standard methods. 

A typical simulation proceeds as follows: (i) Choose the necessary input parameters: fi, T, 
atomic species, trapping potential; (ii) Create an ensemble of realisations, each corresponding 
to a unique set of noise realisations at each time step; (iii) Propagate this set of realisations to 
equilibrium, i.e. when observables (e.g. condensate number) become constant on average; equi- 
librium observables may be extracted by constructing correlation functions from the set of noise 
realisations, e.g. the density is given by (0*0) = Y/iLi ^*^i/N where i denotes a particular noise 
realisation; (iv) Once at equilibrium, dynamical perturbations may also be studied, e.g. topological 



excitations (Section IVB ), or collective modes. 



E, Interpretation: Single Runs and Extracting Coherence Properties 



1. Single Versus Averaged Runs 

By construction, physical properties are meant to be calculated by averaging over different re- 
alisations of the stochastic field O. Nonetheless important information may also be extracted from 
single numerical runs; in this sense, the SGPE ofi'ers a strong analogy to experimental methods, 
in which data is obtained through repeated measurements from several independent experimen- 
tal realisations. The SGPE was first applied in this way to demonstrate that important details of 
the growth and collapse dynamics of 'Li condensates [|27]| contained within single numerical re- 
alisations were lost by averaging over many runs 13, suggesting that single stochastic numerical 
realisations are analogous to independent experimental realisations. The role of information ex- 
tracted from single runs has been strengthened by further analysis, including spontaneous vortex 



formation via the Kibble-Zurek mechanism [(28l|, fluctuating soliton dynamics (Section [IV B|) and 



in situ density fluctuations in atom chip experiments (Section IV A). 



2. A Posteriori Condensate/Quasi-Condensate Extraction 

The noisy wavefunction of the SGPE represents both coherent and incoherent atoms within 
low-energy modes in a unified manner, and a further statistical analysis is required to iden- 
tify the coherent components. The density coherent quasi-condensate may be identified via 
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FIG. 1. Left Plot: Growth of total (solid) and condensate number (dashed) during equilibration: insets 
show corresponding profiles at indicated times, with equilibrium results (at time fgq) shown by the thin 
solid and dotted lines. Right (Top): Condensate number distribution during growth at indicated times. 
(Bottom): Higher temperature data showing SGPE total (thick black noisy hne) and quasi-condensate (thin 
black noisy) density profiles versus modified Popov total (thick grey lines) and quasi-condensate densities 
(dot-dashed grey). The SGPE Penrose-Onsager condensate (shaded region) is compared to the result of 



Eq. ( 16 1 (thin grey line). Figure adapted with permission from S.P. Cockburn et ah, Comparison between 
microscopic methods for finite-temperature Bose gases, Phys. Rev. A, 83, 043619 (2011) Il22l . Copyright 
(201 1) by the American Physical Society. 



«qc(^) = a/2 - g^^\x) n{x) [[T9l |29l |30l, where the second order correlation function g^^\x) = 
(|0(x)|'*)/(|0(x)p)^. The additionally phase coherent fraction is associated to the Penrose-Onsager 
condensate mode [|3TI . identified as the eigenmode corresponding to the largest eigenvalue of the 
system one-body density matrix p{x,x') ~ (0*(jc)0(x')) [E |22l |32l; this may be numerically 
obtained by diagonalising the density matrix — see Fig. [T]for an analysis of a T > Id Bose gas 
based on this prescription. Motivated by Refs. [fT6l[33ll . the Penrose-Onsager condensate of the 
SGPE is found to be well matched within a trapped system by the definition [12211 

n',{x) = g^'\0, X) ^l2-g^^Kx)n{x\ (16) 

where g^^\Q,x) = p(0,x)/ y/n(0)n{x) is the first order normalised correlation function. The de- 
pendence upon ^^^^(0, x) illustrates clearly the additional phase coherence of the Penrose-Onsager 
condensate, relative to the quasi-condensate in which only density fluctuations are suppressed. 



Eq. (16) provides an alternative means of extracting the phase coherent fraction of a trapped gas 



from SGPE simulations, which accurately captures the condensate edge (but breaks down at very 
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small distances from the trap centre); this method is also ideal for distinguishing between 'con- 
densate' and 'quasi-condensate' in atom chip experiments [|22l . 



3. Comparison to the Modified Popov Method 

As an independent validation of the above interpretation, Fig.[T]compares the SGPE result to the 
modified Popov theory of Stoof and co-workers [[T6l[33] - l35l . which accounts for contributions of 
phase fluctuations to all orders. Densities here are obtained within the local density approximation, 
via 



1 



n{x) = n^cix) + y 2j 

k 



1 



with HqcCx) = [yu - Vsxtix) - 2gun'{x)] /gu, where et = [sj + 2gunqcSkV'^ is the Bogoliubov 
dispersion relation, Sk = ff'k^/2m, and V is the system volume. Building on excellent agreement 
between the SGPE and modified Popov theory [[T6l . and the ideas introduced in 133]. the present 
comparison corroborates the above means of extracting the 'true' and quasi-condensate fractions 



of the gas from SGPE data. Thus, Eq. ( [T6| ) is also expected to be a useful tool for analysing 
experimental density profiles. Further comparison between the SGPE and other one-dimensional 
Rose gases theories may be found in [ , 22,| . 



ni. VALIDITY ISSUES 



A. Validity Domain 



Two main assumptions underlie Eq. ( [TT] ): Firstly, high energy thermal atoms within the system 
are treated as being at equilibrium (with their mean field contribution to O currently neglected). 
Thermal cloud dynamics — crucial when the thermal cloud is strongly perturbed — can be in- 
cluded by evolving the distribution functions via the quantum Boltzmann equation [Eq. (|4])] [|5l. 

Second is the so-called 'classical approximation': This terminology stems from the fact that the 
classical Rayleigh- Jeans distribution arises as the leading order term in a small /^(ec -jX) expansion 
of the Bose-Einstein distribution. While it does not constitute an essential ingredient of the theory, 
this approximation is very useful for numerical purposes, as it simplifies the scattering term R{r, t) 



to the form of Eq. ([T0|), thus leading to the SGPE of Eq. (11) that is numerically solved. Although 



this is a well justified approximation for highly occupied (thus low energy) modes, it does lead 
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to an 'ultraviolet catastrophe', which manifests through a dependence of physical observables 
upon the energy cutoff |l6l|20l|; this problem is far more pronounced in spatial dimensions greater 
than one, due to the form of the density of states, and a possible solution is the introduction of 
divergence-cancelling counter-terms [36J. The stochastic field O thus represents not just the phase- 
coherent part (condensate), but all atoms within modes up to an energy cutoff" (in our simulations 



this is typically set by the spatial discretisation — see also Section IV A). 



B. Relevance to Other Theories 

The SGPE is related to a number of theories discussed in this book: The closest link arises to 
the simple growth SPGPE ||2l[37l|38l, which is very similar in nature, despite their rather distinct 
derivations. Various numerical diff"erences arise in practice — most notably the use (or not) of 
a projector to separate low and high-lying modes. While this may be fundamentally important, 
the potential benefits from its use for dynamical predictions when low-lying modes are coupled 
to a static thermal cloud are not universally accepted — see Appendix C2 of Ref. O for a more 
detailed discussion of the links between these two approaches. 

The SGPE is a grand canonical theory, as the exchange of particles and energy between the 
low-energy modes and heat bath is allowed. If, upon reaching equilibrium, fiL^ is set to zero in 



both damping and noise terms (i.e. iR = j] = O 'm Eq. ([TT])), then it reduces to a multimode, finite 
temperature time-dependent GPE with a stochastically- sampled initial state. Such an approach, 
first applied by us to study quasi-condensate growth on an atom chip [18], and subsequently used 
for finite temperature vortex dynamics ||^ , is similar in spirit (but not implementation Q) to the 
stochastic sampling of the Wigner distribution function evolved in truncated Wigner simulations 

Moreover, classical field methods for Bose gases P3l - l48ll . based upon the observation that 
the GPE accurately describes the dynamics of all highly occupied modes, typically start with a 
suitably random, multimode initial condition; this evolves to a classical equilibrium, sampling 
a microcanonical phase- space under ergodic GPE evolution [|46l - |49ll . Since ji and T are input 
parameters for the SGPE, the latter approach may therefore be viewed as a more controlled way of 
generating a finite temperature initial state for classical field simulations [22] . However, one way 
in which the SGPE contrasts to 'conventional' classical field theories is through the generation of 
an ensemble of independent realisations, with physical observables (such as correlation functions) 
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generated by ensemble averaging over many noise realisations (although see, for example, also 
ll50ll ). In classical field theories, such observables are instead typically generated by sampling the 
system at many different times, chosen so as to be sufficiently far apart. 

By construction, the SGPE incorporates fluctuations into the condensate mean-field stochasti- 
cally, whereas these can at most be a posteriori included in theories based on symmetry-breaking. 
Although we have focused here on numerical realisations with a static thermal cloud, the SGPE 
is in general intended to be coupled to a quantum Boltzmann equation for the thermal cloud; in 
this sense the full theory of Stoof may be considered as an (explicitly [/(l)-symmetry-preserving) 
generalisation of the ZNG method [fT4l[T5ll . The main importance of including fluctuations is for 
describing dynamical processes in the region of critical fluctuations (enhanced in low-dimensional 
systems), or for accounting for experimental shot-to-shot variations, whereas a ZNG-type ap- 
proach could only account for averaged properties, albeit doing so very accurately. Finally, set- 



ting 77 ^ but maintaining the dissipative contribution to Eq. (11) leads to the commonly used 
Dissipative GPE (DGPE), with an ab initio expression for the damping rate, as opposed to phe- 
nomenological input. 



IV. APPLICATIONS 

A. Comparison to Quasi- Id Bose Gas Experiments 

Atom chips [l54l facilitate controlled experiments with weakly-interacting, effectively one- 
dimensional Bose gases [l5T| - |53ll55l - l64l , including in situ measurements, thus allowing for preci- 
sion tests against theory. Importantly, fluctuations play a key role over a wide temperature regime 
for such extremely elongated geometries, thus making the SGPE a prime candidate for modelling 
such systems ||23]| . Figure |2] shows a comparison between the SGPE and several, independent 
sets of quasi- Id experimental data. The agreement is excellent when comparing both density 
profiles [Figj2];a)-(b)] and density fluctuations [Figj2]^d)]; the latter requires the binning of raw 
SGPE data [see Figj2];c)] to bins of width set by the resolution in a given experiment, a key step 
required to achieve a consistent analysis of fluctuations. Experimentally, a system is considered 
as one-dimensional if iJ.,kQT fio)^. If the first requirement is not satisfied, the replacement 
^idli^l^ fi-ioA Vl + 4a|(AP - 1] [[65] - [67l can account for the transverse swelling of the gas; in 
addition, ifk^T <jt ha)j_, then we should also account for atoms in transverse excited modes which 
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X [pixels] <N> 

FIG. 2. Quasi-ld SGPE model (solid black line) versus experimental density profiles (circles) from (a) 
Trebbia et al. ISTl and (b) van Amerongen et al. |52]. (c) Main Plot: Single -run density, binned to 
experimental CCD camera resolution for the experiment of Armijo et al. |[53l (grey noisy) vs. average over 
1000 realisations (black), with corresponding raw SGPE single-run data shown in the inset, (d) Density 
flutuations per binned average number predicted by the quasi-ld (black diamonds) and Id (squares) SGPE 
[based on the binning procedure illustrated in (c)] versus experimental data of Armijo et al. (circles). Figure 
adapted with permission from S.P Cockburn et al., Quantitative study of quasi-one-dimensional Bose gas 
experiments via the stochastic Gross— Pitaevskii equation, Phys. Rev. A 84, 023613 (2011) Il23l . Copyright 
(201 1) by the American Physical Society. 

contribute a density = (l//ldB) Z7=iO+l)gi/2(e^~^=^'^''^"^''*''^^^*^^^) [Ell (gi/2(»): polylogarithm 
of order 1 /2). These two amendments yield a quasi-ld SGPE [23J, which is cutoff independent (as 
both below and above cutoff physics is included in the model in an approximate, but self-consistent 
manner), and thus accurately models experiments in the crossover from one to three dimensions 

mm. 

B. Dark Soliton Dynamics in a Quasi-Condensate 

As the equilibrium state of the SGPE agrees well with both experiment ||23|| and alternative the- 
ories in suitable limits [|22|. it constitutes an ideal ab initio approach for finite temperature Bose 
gas dynamics [m EB El [6l (see also [[171 |28l |39l EQl El)- This is particularly true for pertur- 
bations which do not push the thermal cloud far from equilibrium, making the dynamics of a dark 
soliton a perfect candidate. Including fluctuations in the background field leads to 'shot-to-shot' 
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FIG. 3. (a) Histogram of soliton decay times and (b) example SGPE soliton trajectories at T = 175nK, from 
the indicated bins, shown by vertical dot-dashed, solid and dashed black lines; the solid grey trajectory 
corresponds to the DOPE result, (c) Average decay times extracted from soliton decay time histograms 
at several temperatures, the fits to which are shown in the inset. Subplots adapted with permission from 
S.P. Cockburn et al., Matter-Wave Dark Solitons: Stochastic versus Analytical Results, Phys. Rev. Lett. 
104, 174101 (2010) |21j | and S.P. Cockburn et al., Fluctuating and dissipative dynamics of dark solitons in 
quasi-condensates, Phys. Rev. A 84, 043640 (2011) [24]. Copyright (2010) and (2011) by the American 
Physical Society. 



variations in soliton behaviour, as evident from indicative trajectories shown in Fig. [3];b). The 
corresponding histogram of decay times over the ensemble of realisations (Fig. [Sj^a)) is well-fitted 
by a lognormal distribution, displaying an extended tail at long times, indicative of very long-lived 
solitons, relative to the average soliton decay time [l2Tll24ll69l , consistent with experiments 
These distributions are shown for a range of temperatures in Fig. [Sj^c) (inset), with the average 
times found to vary with temperature as (r) ~ T""* (main plot). The distribution of soliton decay 
times obtained via the SGPE, relative to the single mean field DGPE result, indicates that consider- 
ation of many realisations of the stochastic wavefunction allows one to construct a representation 
of the full probability distribution for the gas; moreover, our analysis highlights once more the 
important feature that useful information is also retained within individual SGPE runs, which thus 
bear strong similarities to single-shot experimental realisations [r71 l2Tll23ll24ll28ll69l . 
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